Bistability in superconducting rings containing an inhomogeneous Josephson junction 
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We investigate the magnetic response of a superconducting Nb ring containing a ferromagnetic 
PdNi Josephson junction and a tunnel junction in parallel. A doubling of the switching frequency is 
observed within certain intervals of the external magnetic field. Assuming sinusoidal current-phase 
relations of both junctions our model of a dc-SQUID embedded within a superconducting ring 
explains this feature by a sequence of current reversals in the ferromagnetic section of the junction 
in these field intervals. The switching anomalies are induced by the coupling between the magnetic 
fluxes in the two superconducting loops. 



PACS numbers: 85.25.Dq 

A Superconducting Quantum Interference Device 
(SQUID) consists of a superconducting loop interrupted 
by one (rf SQUID), or two (dc SQUID) Josephson junc- 
tions. The magnetic moment of rf SQUIDs and the 
critical current of dc SQUIDs are periodic functions of 
the magnetic flux $ enclosed by the loop. In most 
cases, the periodicity is given by the magnetic flux quan- 
tum $o = h/2e. This is a consequence of the gauge- 
invariant connection between $ and the phase difference 
ip across the junction(s) via the current-phase relation 
(CPR) I s ((fi). Under certain conditions, the CPR of the 
junction(s) is not a simple sine, but can contain higher 
harmonics: I s (^>) = Y] n Ic" sinrup, where the coupling 

(n) 

coefficients I c quantify the relative strength of pro- 
cesses with a coherent transfer of n Cooper-pairs pQ . In 
some cases, the first order coefficient Ic can even van- 
ish, for instance, in asymmetric 45 degree grain bound- 
ary junctions in d-wave symmetry superconductors [2J[3], 
in out of equilibrium SNS junctions |J and for ballistic 
SFS junctions at the 0-tt transition [5, 6 . The remain- 

(2) 

ing second order coefficient Ic will result in a doubling 
of the frequency in the interference pattern. Although, 
a frequency doubling has been observed experimentally 
[7], its origin is still under debate, since there are also 
dynamic effects in inhomogeneous junctions, which can 
lead to this effect [5]. 

Recently, it was observed that the critical current of 
high temperature superconductor dc-SQUIDs shows a 
<i>o/2 periodicity in certain sections of the interference 
pattern [3] . To explain this observation it was suggested 
that the random faceting of the grain-boundary induces 
a distribution of and n couplings along the junction. 
For some particular values of the applied magnetic field 
the first order coupling of the overall junction is zero, 
allowing second order coupling to be dominant. A sim- 
ilar effect has been also measured in the multi terminal 
transport of a 2-dimensional electron gas connected to a 
superconducting loop [9j. 

Here we report on a different material system, show- 



ing a similar phenomenology and suggest an alternative 
interpretation of this phenomenon. We study the mag- 
netic response of an rf SQUID where the junction is in- 
homogeneous and formed by the parallel connection of a 
conventional (0) Josephson junction and a ferromagnetic 
7r junction. By measuring the total flux in the SQUID 
while increasing the external magnetic field B extl we 
mostly observe a $o-periodic penetration of flux quanta 
into the loop, every time the critical current of the junc- 
tion is exceeded. However, for some values of the ex- 
ternal magnetic field we find a doubling of the switch- 
ing frequency. The proposed model explains this effect 
in terms of a bistability of the supercurrent in the 7r- 
junction for certain values of the applied magnetic field. 
Our double SQUID model successfully explains the mag- 
netic field and temperature dependence of our observa- 
tions, although the CPR of both junctions is assumed to 
be purely sinusoidal. 

We use Nb as the superconductor and dilute PdNi as 
the ferromagnet for our SFS junctions. To pattern the 
loops, we use a robust SiaN^PES mask system for the 
shadow evaporation. The PES (polyether sulfone) forms 
a highly thermostable sacrificial layer [ID]. The 60 nm 
thick Si3N 4 was deposited by plasma enhanced chem- 
ical vapor deposition on top of the PES and provides 
sufficient mechanical stability to resist the large stresses 
created by the Nb film. After patterning the mask by 
electron beam lithography and reactive ion etching using 
CHF 3 , the SiaN4 mask was underetched by an isotropic 
oxygen-plasma. The undercut can have a value of up to 
1 /zm. Evaporation of 40 nm of Nb and 10 nm of PdNi un- 
der different angles provides the desired superconducting 
loops with an integrated SFS planar junction, as illus- 
trated by a scanning electron micrographs of a sample in 
Figs. [T]i and b. In Figs. [IJ; and d the equivalent schemat- 
ics are shown. The thickness of the PdNi film was chosen 
to produce a tt junction close to the 0-7T crossover 
A 10% misalignment of the sample during evaporation of 
the second Nb-layer resulted in an overlap of the two Nb- 
films without PdNi-interlayer, as indicated by the arrow 
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FIG. 1: (a) Scanning electron image of the loop on top of the 
active area of the Hall-sensor. Two leads connected on both 
sides of the junction are used to control the flux through the 
loop, (b) A zoom onto the junction area showing the parallel 
connection of the Nb/Nb-contact (junction A, vertical lines) 
and the Nb/PdNi/Nb- (junction B, horizontal lines), (c) 
Schematic of the device geometry, (d) Schematic of the equiv- 
alent double SQUID: Loops 1 and 2 with self-inductances L\, 
L 2 and mutual inductance L i2 containing junctions A and B. 



in Fig. [T]d. Strong gettering of residual gas by the Nb 
during evaporation of PdNi results in a rather transpar- 
ent tunneling contact between the two Nb-films in this 
area, as a second Josephson junction. 

We have placed the sample on top of the active area of 
a micron-sized Hall-sensor in order to detect the magnetic 
response [T5J [T3] • The Hall-sensor is realized in a semi- 
conductor heterostructure having the electron density of 
2.25 x 10 15 m- 2 and the mobility of 1.13 x 10 6 cm 2 /Vs. 
We achieve a sensitivity of roug hly 500nT/Hz 1/2 to 
SOnT/Hz 1 / 2 depending on the sensor current. For our 
loop dimensions of approximately 7.6 /im x 8.5 /xm, the 
magnetic flux quantum $o = 2.067 x lCP 15 Vs corre- 
sponds to a magnetic field of about 31 /iT. 

Upon sweeping the external magnetic field, circulating 
supercurrents in the loops are generated. The resulting 
flux in the loop, $i, induces a flux-periodic contribution 
to the Hall voltage across the Hall sensor, while the con- 
tribution of the external flux $ e i is simply subtracted. A 
typical trace of the induced flux Li/i = <f>i — <f> e i vs. <i> e is 
depicted for two different temperatures in Fig. [2^. Only 
the flux through the ring is measured due to negligible 
contribution of the almost orthogonally tilted micro loop, 
see Fig. [TJ;. The magnetic response of the ring is strongly 
hysteretic due to the large L/ c -product, the latter is char- 
acterized by the parameter /3li = 2ttLiI c i/^q 3> 1 [T3] . 
The ring inductance is determined from the estimated 
filling factor of the magnetometer. The critical current is 
determined from the vertical size of the hysteresis loops 
in Fig. |2^. 

The signal shows additional substructures in the 
switching pattern in certain intervals of the external mag- 
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FIG. 2: (a) Full hysteresis cycle for two different T = 1.5 K 
(outer cycle) and 3.0 K (inner cycle), corresponding to the 
critical currents I c \ ~ 980 and 690 \xA respectively. The 
arrows indicate positions of substructures with double fre- 
quency. With increasing T the substructures are shifted to 
higher fields. The ring inductance L\ ~ 26 pH is determined 
from the estimated filling factor of the magnetometer. The 
irregular drifts and jumps are artifacts of the Hall cross, (b) 
Magnified top substructure on outer cycle is shown as an ex- 
ample (All other substructures are similar). The absolute 
value of the induced flux $i — $ e i is determined from the 
jump height with an uncertainty of 0.1$o- 



netic field, as indicated by the arrows in Fig. pp. This 
effect has been seen in three samples with similar inho- 
mogeneous junction geometry. A zoom onto the top sub- 
structure in the outer cycle is shown in Fig. [2}). All other 
substructures are similar. It can be seen that the regu- 
lar $o-periodic switching pattern is interrupted by addi- 
tional peaks, which gain height, until they take over. The 
field intervals displaying the substructure in the switch- 
ing behavior shift towards higher field when the tem- 
perature is increased. The substructure looks similar to 
the predicted frequency doubling effect, expected for a 
dominant second harmonic contribution to the CPR |6J. 
However, as we show below, the observed substructure 
can be explained in terms of two coupled loops, even for 
a sinusoidal CPR of both and 7r junctions. 

In the following model treatment, we approximate the 
extended inhomogeneous junction by a parallel connec- 
tion of two short junctions A and B, the latter being in 
the 7r-state. We assume a sinusoidal current-phase rela- 
tions Ia(b) — IcA(B) sm <PA(B)- These two junctions form 
a small dc SQUID with an inductance L2, which inter- 
rupts the large loop with inductance L\. As sketched in 
Fig. [TJl, our model system then consists of a rf SQUID, 
with an embedded dc SQUID as weak link. The super- 
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currents circulating in the two SQUID loops are coupled 
by the Kirchhoff laws and their mutual inductance Li 2 . 
The free energy of the circuit is given by 

W = E A {^ A ) + E B {i PB )+ l -L 1 I 2 1 + l -L 2 I 2 2 +L 12 I 1 I 2 (1) 

where tpA and ifs are the macroscopic phase differences 
across the junctions A and B, I\ = I a + Ib is the cur- 
rent circulating in the large loop, I 2 = I a is the cur- 
rent circulating in the small loop, where 1a and Ib are 
the currents through junctions A and B, respectively. 
The first two terms in Eq.([T|) are the Josephson ener- 
gies Ei(<pi) = ($ /27r)|I ci |(l T cos<^), (i = A, B) for Gl- 
and 7r-junctions, respectively. The three remaining terms 
represent the magneto-static energy [14j [15] . The mag- 
netic fluxes $1 and <& 2 through the loops 1 and 2 are 
given by 

$i = $ e i-Li/i-Li 2 / 2 , $ 2 = <S>e2-L 12 I 1 -L 2 I 2 . (2) 

Here, $ e i and <fr e2 are the corresponding fluxes of the 
external magnetic field. The total fluxes are related to 
the phase differences 

*1 = ^^B, $2 = ^(<PA - fs)- (3) 

Finally, Eq.([T]) can be rewritten explicitly as a function 
of i»i and $2 in the form 



W = E a (va) + Eb(vb) 
L\L 2 — L\ 2 { 2 

- Li 2 ($l - $el)($2 - *e2)}- 



Li 



(*2 



1>e2) 2 

(4) 



For given external fluxes <& e i and <i> e 2, the local min- 
ima of the free energy W with respect to <&! and <&2 are 
calculated numerically and plotted in Fig. [3] Like for 
usual rf SQUIDs in the hysteretic regime, the circulat- 
ing current in the large loop is a multivalued function of 
external flux if /3li > 1- When the external magnetic 
field is swept, e.g., in negative direction, the circulating 
supercurrent follows the lines of local energy minima, as 
indicated by the red lines in Fig. [3] At the upper end 
of each line the state becomes unstable and the system 
switches into the nearest available state with lower free 
energy, as reflected by the sharp drops in <i>i. The en- 
velope of these jumps is periodic, due to a modulation 
of the maximum I\ by the integrated dc SQUID. This 
comes from the coupling between <!>i and $2, implying 
a circulating current also in the small loop (e.g. with 
Ib > 0), which contributes to the free energy. Around 
^ei/^o ~ 5, states with opposite current Is in the small 
junction become stable (see the short green lines in the 
upper inset in Fig. [31). The stability region for this set 
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FIG. 3: The map of local energy minima of the double 
SQUID. Red and green lines correspond to positive and neg- 
ative sign of the current Ib- Top and bottom sawtooth lines 
represent the largest hysteresis in $1 — $ e i. Parameters are 
the same as in Fig. [4] (outer cycle). Upper inset: Nucleation 
of states with opposite (green). Lower inset: Alternation 
of states (red and green lines) with positive and negative sign 
of Ib in the region with frequency doubling. 



of states grows until they eventually become more sta- 
ble than those with the original (positive) orientation of 
Ib- This is precisely the region, where these states be- 
come involved in the switching process. In this region, 
states with positive and negative orientation of the cur- 
rent in the small loop alternate, resulting in a doubling of 
the switching frequency (lower inset in Fig. [3]). When the 
external flux is decreased further, the states with the pos- 
itive orientation of Ib become energetically unfavorable 
and their region of stability shrinks, until the switching 
processes entirely involve states with negative Ib (green 
lines). In this way, the doubling of the switching fre- 
quency is traced back not to a period doubling in the 
CPR, but to the presence of the two-fold orientation of 
the current in the 7r-section of the junction. 

To compare the numerical results with our experi- 
ment we take junction A in the state (tunnel junction, 
IcA > 0), junction B in the it state (SFS, I c b < 0), and 
Ica/IcB = -10, hi - 90, 0L2 = 2irL 2 I cB /<f> = -1.5, 
27tLi2/ c a/*o = -9.7 and $ e2 /$ e i = -0.005. The crit- 
ical current I c a and inductance of the ring L± are de- 
termined directly from the experimental hysteresis loop 
in Fig. [2J1, while I c b is estimated from previous measure- 
ments on similar Nb-PdNi junctions |12j . The dc SQUID 
inductance L 2 and mutual inductance L12 affect mainly 
size and position of the switching anomaly on the hys- 
teresis loop. For reasonable values of chosen parameters 
a good agreement is achieved between the experimental 
data for T = 1.5K and numerical simulations (see Figs. 
[2] and [4]). Note that the ratio $ e2 / < I ) e i of external fluxes 
differs from the loops area ratio and has (for this sam- 
ple) negative sign due to the almost orthogonal tilt of 
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FIG. 4: (a) Calculated largest hysteretic cycles fitted to the 
experimental data shown in Fig [2] with the parameters given 
in the text. The inner cycle corresponds to a reduction of I c a 
and I c b by 25%, corresponding to the higher temperature. 
The dotted lines correspond to the reversal of the sweep direc- 
tion in the experiment. Inset: Shift of the switching anomaly 
with the measured total critical current Ic(T) at tempera- 
tures 1.5, 2, 2.5 and 3 K (dots). The solid line represents 
the theoretical prediction with the same set of parameters 
(b) Zoom on the upper substructure in the outer cycle, (c) 
Corresponding induced flux through the small loop. 



the small loop with respect to the sample. Its value de- 
termines the number of the observed double switching 
events and it can be readily determined from the data. 

If we assume that the measured reduction of I c a to 
« 75% at T = 3 K (see Fig. [4} is similar also in \I c b\, 
we can well reproduce the observed shift of the switch- 
ing anomaly towards higher external flux (inner trace in 
Fig. |4^i). The dots in the inset in Fig. [4^ show the shift of 
the switching anomaly vs. the critical current Ic{T) for 
different temperatures. The solid line denotes the pre- 
diction of our theory using the same set of parameters. 

Sequences of the frequency doubling in the flux mod- 
ulation periodically occur in certain intervals of the ex- 
ternal magnetic field, if \J3ia\ > 1 of the smaller loop 
is sufficiently large. A zoom onto the substructure is 
shown in Fig. [4)3. Substructures in the switching pat- 
tern $1 — <& e i occur with the period |$ei/3>e2|- The cor- 
responding flux modulation in the small loop is shown in 
Fig. [2J;. When the external field is varied, small periodic 
perturbations induced by the large loop alternate period- 



ically the energies of the two opposite current directions 
in the small loop, thus forming the bistable region with 
large flux oscillations. The bistable region is placed sym- 
metrically around the external flux value corresponding 
approximately to integer number of $0 m the small loop 
and equal energies for opposite currents in junction B. 

Additional calculations show that in the ground state 
the first substructure occurs in the low external field 
corresponding to <I> e 2 ~ [16] . When both junctions 
arc in the state the first substructure is located at 
<J> e2 ~ < &o/2, which corresponds in our case to a very 
large «3> e i. However, positions of substructures on the 
hysteretic loop are strongly shifted from the ground state 
values, and depend on the mutual inductance. Numeri- 
cal calculations show that practically the same hysteretic 
behavior shown in Fig. |4]can be obtained for junctions A 
and B both in the state with similar parameters of the 
double SQUID [TB]. 

To conclude, we studied an rf SQUID containing an 
inhomogeneous Josephson junction as a weak link. We 
have found experimentally a doubling of the switching 
frequency in certain ranges of magnetic flux. The in- 
homogeneous junction can be modelled as a small dc 
SQUID with and 7r Josephson junctions. This model 
explains the observed switching anomaly by a bistable 
switching of the orientation of the current in the weaker 
section of the junction. The suggested mechanism is ef- 
fective, independent of the shape of the CPR in both 
junctions, and may also be relevant for similar observa- 
tions in other systems. 
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